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[l] We investigate the transport of ions in the low-altitude 
magnetosphere of Mercury. We show that, because of small 
spatial scales, the centrifugal effect due to curvature of the 
E x B drift paths can lead to significant particle energiza- 
tion in the parallel direction. We demonstrate that because 
of this effect, ions with initial speed smaller than the escape 
speed such as those produced via thermal desorption can 
overcome gravity and escape into the magnetosphere. The 
escape route of this low-energy exosphere originating 
material is largely controlled by the magneto spheric con- 
vection rate. This escape route spreads over a narrower 
range of altitudes when the convection rate increases. Bulk 
transport of low-energy planetary material thus occurs 
within a limited region of space once moderate magneto- 
spheric convection is established. These results suggest that, 
via release of material otherwise gravitationally trapped, the 
E x B related centrifugal acceleration is an important 
mechanism for the net supply of plasma to the magneto- 
sphere of Mercury. Citation: Delcourt, D. C., K. Seki, N. 
Terada, and T. E. Moore (2012), Centrifugally stimulated exo- 
spheric ion escape at Mercury, Geophys. Res. Lett., 39, L22105, 
doi: 10.1 029/201 2GL054085 . 

1. Introduction 

[ 2 ] Following the pioneering observations of Mariner-10 
in 1974-75, in-situ measurements from the MESSENGER 
spacecraft confirm the existence of an intrinsic magnetic 
field at Mercury [e.g., Anderson et al., 2011] ; hence, the 
formation of a miniature magnetosphere with spatial and 
temporal scales much smaller than those at Earth (by about a 
factor 8 and a factor 40, respectively). The typical standoff 
distance of Mercury’s magnetopause is expected to be of 
the order of 1.5-2 R M [e.g., Slavin et al., 2010a, 2010b]. 
Material of planetary origin can be released into the exo- 
sphere of Mercury via a variety of processes such as thermal 
desorption, electron-stimulated desorption [e.g., McLain 
et al., 2011], photon-stimulated desorption, solar wind 
sputtering or micro-meteoritic bombardment [e.g., Leblanc 
and Johnson, 2003], Once ionized by solar UVs, these 
newly formed ions of planetary origin are subjected to a 
large-scale transport and contribute to the different regions of 
Mercury’s magnetosphere, as found from MESSENGER 
observations and simulation modeling [e.g., Zurbuchen et al, 
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2008, 2011; Delcourt et al., 2003; Paral et al, 2010], The 
purpose of this study is to examine the impact of centrifugal 
effects in the miniature environment of Mercury on the fate of 
exospheric material. These centrifugal effects will first be 
described in section 2. In section 3, we will examine their 
consequences on the transport of low-energy ions. 


2. Centrifugal Acceleration 

[ 3 ] In the adiabatic (guiding center) approximation, the 
parallel equation of motion in the absence of parallel electric 
field writes [see, e.g., Northrop, 1963, equation (1.20)]: 
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Here, gy is the gravitational acceleration in the parallel 
direction, \ E is the E x B drift velocity, V\\, the particle 
parallel speed, b, a unit vector in the B direction, .v. the 
curvilinear abscissa along the field line, and m and p, the 
particle mass and magnetic moment, respectively. Equation (1) 
displays two distinct curvature related acceleration terms. 
The third term on the right hand side of (1) is due to curvature 
of the magnetic field lines and can lead to quite substantial 
acceleration upon crossing of a field reversal. The fourth 
term on the right hand side of (1) is due to curvature of the 
E x B drift paths, a centrifugal acceleration that can be 
viewed as that affecting a bead (the charged particle) on a 
rotating rod (the convecting magnetic field line) [e.g., 
Horwitz et al, 1994], In the Earth’s magnetosphere, it was 
shown that this latter E x B related centrifugal effect can 
lead to upward ion energization of a few e Vs or a few tens of 
eVs at low altitudes [e.g., Cladis, 1986] or to focusing 
(referred to as “centrifugal trapping”) of slow plasma sheet 
particles into the equatorial region [e.g., Delcourt et al., 
1990], Because of small spatial scales, it was shown by 
Delcourt et al. [2002] that this E x B related centrifugal 
acceleration is more pronounced at Mercury than at Earth, 
possibly yielding energization of planetary material up to 
several hundreds of eV within minutes when ions are trans- 
ported into the magnetospheric lobes over the polar cap [see, 
e.g., Delcourt et al., 2002, Figure 1]. Assuming a steady 
magnetic field configuration, this acceleration occurs in an 
adiabatic (magnetic moment conserving) manner since, in 
this region of space, the magnetic field does not significantly 
change within a cyclotron turn. As shown in Delcourt et al. 
[1990, section 3.2], an equation similar to (1) can be 
obtained for the full particle with the exception of the mirror 
force term (last term in (1)). The centrifugal acceleration 
described above thus does not depend upon the adiabatic or 
nonadiabatic character of the motion. 

[ 4 ] The centrifugal acceleration due to curvature of the 
E x B drift paths is the one of interest in this study since it 
opposes the downward oriented gravitational acceleration 
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Figure 1. (left) Na + trajectories in the X-Z plane, (right) From top to bottom, Na + altitude, latitude, parallel speed, E x B 
drit speed, and parallel speed-to-E x B drift speed ratio as a function of path length. The test ions are launched from the 
planet surface at 65° latitude in the day side sector, with 0.01 eV or 0.1 eV energy, 170° pitch angle and arbitrary gyrophase. 
Blue, green and orange-red colors correspond to cross-polar cap potential drops of 2 kV, 5 kV, and 10 kV, respectively. 


(g|| in (1)). It can thus affect the net upward transport of ions 
depending upon magnitude of the magnetospheric convec- 
tion. To quantitatively investigate this effect with the help of 
numerical simulations, we adopt as in previous studies [e.g., 
Delcourt et al., 2003] the model of Luhmann and Friesen 
[1979] that consists of the superposition of a dipole field 
with a Harris sheet. As suggested by recent MESSENGER 
measurements [Anderson et al., 2011], we also take into 
account a northward offset of 484 km for the magnetic field, 
and we use the following model parameters : B T = 70 nT, 
L = 0.5 R m , and B E q = 195 nT (B r being the asymptotic tail 
field, L, the half-thickness of the current sheet and B E q, the 
equatorial surface field). As for the large-scale convection 
electric field, we use in first approximation a two-cell Volland- 
Stem pattern [e.g., Volland, 1978] that leads to antisunward 
transport in the polar cap region and sunward transport 
below. Distinct cross-polar cap potential drops are considered 
as described hereinafter. As for the test particle trajectories, 
they were computed using the full equation of motion, inte- 
grated by means of a fourth order Runge-Kutta technique. 

[5] Figure 1 shows an example of the ion trajectories 
obtained. We consider here the sodium population that is 
released into Mercury’s exosphere at very low energies (less 
than 0.1 eV) as expected for instance from thermal desorp- 
tion [see, e.g., Leblanc and Johnson, 2003, Figure 3], Once 
ionized by solar UVs, the newly produced Na + ions are 
entrained by the convecting magnetic field lines and circu- 
late inside the magnetosphere. The trajectories coded in blue 
in Figure 1 (left) show the results obtained for 0.01 eV and 
0.1 eV Na + originating from the planet vicinity (5 km alti- 
tude) in the high-latitude dayside sector. In this case, a weak 
convection electric field is considered, with a cross-polar 
potential drop of 2 kV. It can be seen in this panel that, due 
to initial speeds smaller than the escape speed, the ions are 


trapped by the gravitational field, viz., they travel upward 
until some maximum altitude at which point they are turned 
back toward the planet and subsequently precipitate onto the 
surface (note here that during this transport, the test Na + with 
0.1 eV initial energy is convected into the nightside sector 
over the polar cap). This gravitationally bound behavior can 
be better appreciated in Figure 1 (right) that show the ion 
position and velocity as a function of path length. It can be 
seen that the Na + with 0.01 eV and 0.1 eV energy reach 
maximum altitudes of ^150 km and ~700 km, respectively. 
As expected, their parallel velocity (third panel from top) 
first decreases, then increases upon return toward the planet. 

[6] The trajectories coded in green and orange-red in 
Figure 1 were obtained considering larger cross-polar cap 
potential drops (viz., 5 kV and 10 kV, respectively). In both 
cases, a behavior radically different than trapping at low 
altitudes (blue trajectories) can be seen. That is, regardless of 
initial energy, the test Na + are now transported over the polar 
cap into the nightside sector where they eventually interact 
with the magnetotail current sheet. Not surprisingly, during 
the lobe trajectory sequence, the adiabaticity parameter k 
that is defined as the square root of the minimum curvature 
radius-to-maximum Larmor radius ratio [e.g., Biichner and 
Zelenyi, 1989] is larger than 10 (its initial value 
being > 90) and, in this region of space, the ion motion is 
thus adiabatic (magnetic moment conserving). In contrast, 
during interaction with the magnetotail current sheet, the 
ions are subjected to large magnetic moment change and 
nonadiabatic energization (of several keVs). These ions 
substantially contribute to the plasma sheet populations prior 
to traveling back toward the planet and precipitating in the 
dusk sector. Figure 1 (left) shows that the green and orange- 
red trajectories intercept the pole at altitudes significantly 
higher than the blue ones. As demonstrated by Delcourt et al. 
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Figure 2. Color-coded altitude (a-d) and parallel energy (e-h) of Na + ions at entry into the nightside sector as a function of 
initial energy and altitude. The test ions are launched from the dayside sector at 65° latitude with 170° pitch angle and arbi- 
trary gyrophase. Different cross-polar cap potential drop are considered from left to right: 1 kV, 2 kV, 5 kV, and 10 kV. The 
grey shaded area corresponds to ions that do not cross the pole. 


[2002], the reason for this is that their parallel speed rapidly 
increases after ejection into the magnetosphere under the 
effect of a large E x B related centrifugal effect. Upon 
crossing of the pole, their altitudes are of the order of 0.5 R M . 
Clearly, the low-energy Na + now rapidly overcome the 
gravitational field and can travel up to large distances in the 
Hermean magnetosphere. 

[ 7 ] It is also of particular interest to note in Figure 1 that, 
although two distinct cross-polar cap potential drops are 
considered for the trajectories coded in green and orange- 
red, the altitudes at pole crossing as well as the trajectory 
apex in the nightside sector are fairly similar. Looking at 
Figure 1 (right), it can be seen that the E x B drift speed 
along the orange-red trajectories is twice larger than that 
along the green ones, which should thus lead to a different 
drift path. However, the parallel speed obtained from E x B 
related centrifugal acceleration also increases in the same 
proportion (Figure 1, bottom right) so that the net ion drift 
path is not much affected. In other words, Figure 1 demon- 
strates that the E x B related centrifugal acceleration can be 
significant enough for ions to overcome gravitational trap- 
ping and to escape into the Hermean magnetosphere. 
Moreover, their escape route appears to depend little upon 
E x B magnitude since the ion parallel speed increases in 
proportion to it. Given the large amount of planetary mate- 
rial that is expelled at low energies from the planet surface, 
this centrifugal escape mechanism may play a significant 
role in the net supply of plasma to the magnetosphere. 

3. Centrifugal Escape 

[8] To better characterize the above centrifugal escape 
mechanism, we performed systematic trajectory computa- 
tions for Na + with different initial energies (0.001 eV- 
10 eV) and initial altitudes (below 200 km). The results of 


these computations are presented in Figure 2. Here, to 
explore the large variety of ion behaviors (from gravita- 
tionally bound, as with the blue trajectories in Figure 1, to 
mostly bunched together during escape, as with the orange- 
red trajectories), we consider a wide range of cross-polar cap 
potential drops, from extremely weak (1 kV) to fairly large 
values (10 kV). This leads to a dawn-to-dusk convection 
electric field in the equatorial magnetotail that varies 
between ~0.3 mV. in 1 and ~3.0 mV.m '. The test Na + in 
Figure 2 are launched from 65° latitude in the dayside sector 
on the noon meridian. For those ions that escape into the 
magnetosphere, Figures 2 (top) and 2 (bottom) show the 
color-coded altitude and parallel energy at pole crossing or, 
equivalently, upon entry into the nightside sector. It can be 
seen in Figures 2a and 2e that, for very weak magneto- 
spheric convection, most of the test Na + considered return to 
the planet before crossing the pole. As expected, only those 
ions with initial speed comparable to or larger than the 
escape speed (corresponding to an initial energy of ~2 eV as 
indicated by the solid black line) gain access to the nightside 
sector. 

[ 9 ] On the other hand, if a larger convection rate is con- 
sidered (from left to right in Figure 2), the amount of ions 
that enter the nightside magnetosphere significantly increa- 
ses. As a matter of fact, for potential drops of 5 kV and 
above, all test ions in the initial altitude and energy ranges 
considered escape into the magnetosphere (Figure 2). In 
particular, the comparison of Figures 2g and 2h reveals that 
the Na + parallel energy at pole crossing increases when the 
convection rate increases, as expected from a more pro- 
nounced E x B related centrifugal acceleration. However, as 
discussed in Figure 1, this enhanced parallel acceleration 
does not translate as widely different altitudes upon crossing 
of the pole. As a matter of fact, comparison of Figures 2c 
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Figure 3. Envelope of Na + altitudes at entry into the nightside sector as a function of initial altitude and for different (color- 
coded) cross-polar cap potential drops. Ions with initial energy smaller than (left) 0.1 eV and (right) 1 eV. 


and 2d shows that, as the convection electric field increases, 
the range of altitudes at pole crossing gradually decreases. 
This suggests that, as a result of E x B related centrifugal 
focusing, there exists some limiting trajectory toward which 
the drift paths of all low-energy ions converge. 

[to] This effect can be better appreciated in Figure 3 that 
shows the envelope of Na + altitudes at pole crossing as a 
function of initial altitude and for all energies below 0.1 eV 
(Figure 3, left) and below 1 eV (Figure 3, right). The dif- 
ferent colors in this figure correspond to distinct convection 
rates. It can be seen in Figure 3 that, as the cross-polar cap 
potential drop increases, the range of altitudes at which the 
low-energy ions enter the nightside sector significantly 
shrinks until some narrow channel is obtained (black lines in 
Figure 3). In other words, the pattern in black in Figure 3 
delineates some asymptotic regime in which all low-energy 
ions released at low altitudes overcome gravity and travel 
into the nightside magnetosphere. 

[n] The pattern achieved in Figure 3 was obtained for an 
initial latitude of 65°. For different initial latitudes, patterns 
similar to that in Figure 3 are obtained but at different alti- 
tudes. Accordingly, it can be speculated that the outermost 
altitudes at which low-energy ions originating from the 
dayside enter the nightside sector is associated with a notable 
change of density, a feature reminiscent of the “geopause” 
suggested by Moore and Delcourt [1995] to delineate the 
domain of planetary plasma contribution in the terrestrial 
magnetosphere. This can be better appreciated in Figure 4 


that presents the results of systematic trajectory calcula- 
tions, considering different initial latitudes in the dayside 
sector. Here, a number of test Na + were launched in the 
0-200 km altitudinal range (by steps of 1 km) above 50° 
latitude (by steps of 1°) in the dayside sector, with distinct 
pitch angles and energies (less than 0.1 eV). Once the tra- 
jectory calculations were performed, the density and average 
energy in a given altitudinal bin were calculated by summa- 
tion over the various Na occurring in that bin, in a like 
manner to Delcourt et al. [2003] (see equations (2)— (3) of 
that paper). In Figure 4 (left) that shows the ion altitude at 
pole crossing versus initial latitude, it is clearly apparent that 
the larger the convection rate, the smaller the altitudinal 
spread for given initial latitude. In other words, enhanced 
E x B related centrifugal acceleration yields enhanced 
structuring of the transpolar plasma flow. In terms of density, 
it is apparent from Figure 4 (middle) that this translates as a 
more abrupt variation at high altitude. That is, whereas the 
ion density gradually decreases with altitude for weak con- 
vection (blue profile), a more “step-like” behavior is obtained 
for enhanced convection (red profile). Finally, as expected, 
Figure 4 (right) displays significantly larger total energy 
(solid lines) and parallel energy (dotted lines) at entry into the 
nightside sector when the convection rate increases. 

[ 12 ] It was mentioned above that the MESSENGER 
measurements provide evidence of a northward offset of 
484 km of the magnetic field [Anderson et al., 2011]. This 
offset that was taken into account in the present modeling 



Figure 4. (left) Initial latitude, (middle) density and (right) energy of Na + ions versus altitude at pole crossing. Blue and red 
profiles correspond to two different cross-polar cap potential drops (2 kV and 10 kV, respectively). The dotted lines in 
Figure (right) show the Na + parallel energy. 
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(see section 2), leads to a stronger surface magnetic field at 
high northern latitudes as compared to southern ones. With a 
larger polar cap in the southern hemisphere as compared to 
the northern one, this offset may affect the ejection of exo- 
spheric material in each hemisphere. It should also lead to 
some asymmetry in the ion flows in northern and southern 
lobes. As for the centrifugal escape mechanism, it depends 
upon the magnitude of both convection electric field and 
curvature of the E x B drift paths. This mechanism should 
thus operate in both hemispheres, albeit with different effi- 
ciencies (i.e., different limiting energies for bulk escape) at 
northern and southern latitudes. This efficiency variation is 
clearly reflected in Figure 2, with a change in the amount of 
escaping ions from ~35% up to 100% of the total population 
considered as the convection rate increases. Centrifugally 
stimulated escape may explain the presence of Ca + ions in 
the near-Mercury magnetotail during M3 [Vervack et al, 
2010], a MESSENGER flyby that was characterized by 
intense magnetospheric convection. A thorough assessment 
of this centrifugal mechanism efficiency awaits additional 
in-situ measurements for a variety of magnetospheric con- 
ditions. The MESSENGER measurements also provide 
numerous observations of Flux Transfer Events (FTE) at 
Mercury [e.g., Slavin et al., 2010a, 2010b]. At Earth, such 
localized reconnection phenomena go together with narrow 
channels of enhanced convection at low altitudes [e.g., 
Pinnock et al., 1993]. If this is also the case at Mercury, in 
view of the present results, such FTEs may lead to E x B 
related centrifugal escape of otherwise trapped particles, 
hence ripping off low-energy exospheric material and 
releasing it into interplanetary space. 

4. Conclusion 

[ 13 ] Because of the small spatial scales of Mercury’s 
magnetosphere, the centrifugal acceleration due to curvature 
of the E x B drift paths can lead to significant ion energi- 
zation in the parallel direction. As a result, large amounts of 
low-energy material in the immediate vicinity of the planet 
can overcome gravitational trapping and escape into the 
magnetosphere. The larger the convection rate, the more 
efficient the centrifugal focusing of these escaping ions 
during transport from dayside to nightside, the larger the 
limiting energy below which bulk transport is obtained. This 
E x B related centrifugal mechanism may play a prominent 
role in the net plasma supply to the Flermean magnetosphere 
and/or in the escape of planetary material into interplanetary 
space. 


[ 14 ] Acknowledgments. The Editor thanks two anonymous 

reviewers for their assistance in evaluating this paper. 
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